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Abstract 

We show that a sizable running spectral index suggested by the recent SPT data can be 
explained in the axionic curvaton model with a potential that consists of two sinusoidal contri- 
butions of different height and period. We find that the running spectral index is generically 
given by dus/dhi k ^ (yig — 1), where is the e- folds during one period of modulations. In 
the string axiverse, axions naturally acquire a mass from multiple contributions, and one of the 
axions may be responsible for the density perturbations with a sizable running spectral index via 
the curvaton mechanism. We note that the axionic curvaton model with modulations can also 
accommodate the red-tilted spectrum with a negligible running, without relying on large-field 
inflation. 
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I. INTRODUCTION 



The origin of density perturbations is one of the central issues in cosmology. While 
the quantum fluctuation of the inflaton is the simplest possibility, there may be many 
light scalars in nature, which acquire quantum fluctuations extending beyond the horizon 
during inflation. If so, some of them may give a significant contribution to the observed 
density perturbations via the curvaton jl- J] (or its variant, e.g. modulated reheating js], 
6|) mechanism. Those models can be constrained or preferred by observations, especially 
if one finds an extension(s) of the LCDM cosmology such as non-Gaussianity, running 
spectral index, dark radiation, isocurvature perturbations, and so on. 

Interestingly, the recent SPT data gives preference to a large negative running spectral 
index [3]"^ 

4^ = -0.024 ±0.011 (1) 
dink 

for the CMB (SPT+WMAP7) data alone. The significance increases further if the BAO 
and Ho data are included. In a single-field slow-roll inflation models with a featureless 
potential, the running spectral index is second order in the slow-roll parameters, and 
therefore it is of order (n^ — 1)^ ~ 10~^. Thus, it is a challenge to explain the observed 



running, if taken at face value. There have been various proposals. See e.g. Refs. |9Nl7|. 



In this letter we consider a curvaton scenario [l-4| in order to explain the observed 
running spectral index. We consider an axionic curvaton model in which a pseudo Nambu- 
Goldstone boson plays a role of curvaton. In order to explain the running, we assume 
that the curvaton potential consists of two sinusoidal contributions of different height and 
period. It is possible to generate such potential, if the axion mass receives contributions 



from several instantons, which may be ubiquitous in the string axiverse 18|]. As we shall 



see later, the running spectral index a = dus/dlnk is generically given by 

a (us - 1), (2) 

AN ^ ' ^ ^ 



^ On the other hand, the ACT data [8| does not give preference for any extensions of the LCDM model. 
The SPT and ACT observations look at different patches of the sky, and hopefully the tension will be 
resolved by the Planck or future observations. 



2 



where AA^ is the e-folding number during one period of the modulations. Thus, the 
observed spectral index and its running can be naturally explained if AA^ = (9(10). 

The rest of this letter is organized as follows. In Sec. [TTl we summarize the latest 
observational results and the prediction of the spectral index and its running in the 
curvaton scenario. We present our model in Sec. Illll and show that the observed running 
can be explained in the axionic curvaton model with repeated modulations. The last 
section is devoted for discussion and conclusions. 



II. SPECTRAL INDEX AND ITS RUNNING 



In this section we summarize the latest SPT result and the expression of the spectral 
index Ug and its running a in terms of the curvaton potential V{a)^ which will provide us 
with implications for the curvaton model building. 

The scale dependence of primordial curvature perturbations can be parametrized by 
the spectral index and its running a defined by [3] 

- (3) 

where A^(/c) represents a power spectrum of primordial curvature perturbations TZk-, and 
A; is a comoving wave number. For a constant running, the power spectrum is given in 
the following form, 

A^(A:) = A^(A:o) ^ ■ (5) 



The latest SPT data combined with the WMAP 7yr data give preference to the non-zero 
running at more than 2a [7| 

a = -0.024 ±0.011 (6) 

with ns{k^^^) ^ 0.96 at the SPT pivot scale k^^'^ = 0.025 Mpc"^ If we extrapolate 
the above result to the WMAP pivot scale k^^^^ = 0.002 Mpc~^, the spectral index is 
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roughly estimated to be ns{k^^^^) ~ 1.02 > 1. Thus, the spectral index changes from 
blue-tilted to red-tilted as one goes from k^^^^ to ^q^""". The e-folding number between 
the two pivot scales is ANp-„ot = ln(A;SPT/fcWMAP) ^ 2.5. 

In the curvaton scenario, the spectral index and its running can be expressed in terms 
of the curvaton potential V{(j) and the inflation scale |2oj |: 

-.-1 + 2- + ^, (7) 

~ HI Ht 3 HI HI 9 Ht ' ^ ^ 

where a* and H^ are the curvaton field value and the Hubble parameter when the cosmo- 
logical scale k exited the horizon during inflation. The prime denotes the derivative with 
respect to a. 

Let us make a simplifying assumption that the Hubble parameter hardly evolves during 
inflation, which is justified except for the large-field inflation.^ Then, the above expres- 
sions can be simplified as 

^1 + ^^, (9) 

9 Ht ^ ^ 

In order for Ug to change from > 1 to < 1, therefore, the curvature of the potential 
must change its sign from positive to negative as the curvaton evolves. This requires 
some structure in the curvaton potential. However, if such structure is present only in a 
limited field range, we would encounter another problem: why cosmological scales exited 
the horizon just when the curvaton passed the region where the structure exists. This 
fine-tuning problem can be avoided if the potential has a structure everywhere in the 
potential. Indeed, Kobayashi and the present author applied this idea to the large-field 
inflation by adding repeated modulations to the inflaton potential, and showed that the 
large negative running can be realized without affecting the overall inflaton dynamics |l7| . 



We note that the following argument can be applied to large-field inflation in a straightforward way. 
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This should be contrasted to the resuh of Ref. 16j that the total e-folds of inflation does 
not exceed about 30 in the presence of large negative running. In the same spirit of 
Ref. 171], we consider the curvaton potential with modulations in the next section. 



III. AXIONIC CURVATON 

We consider an axionic curvaton model, in which the curvaton is a pseudo-Nambu- 
Goldstone boson. The shift symmetry is assumed to be explicitly broken by e.g. non- 
perturbative effects as in the QCD axion, generating the following potential: 

,4 ' ^ , a 



V,{a) = A4(^l-cos(^-JJ. (11) 
In the following we will focus on — tt < a / f\ < it without loss of generality. The axionic 



curvaton has been studied in e.g. Refs. 
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-|24|. The model has an advantage that 



the curvaton mass is protected by the (approximate) shift symmetry which enables the 
curvaton to acquire quantum fluctuations during inflation.^ 

The curvaton cr acquires quantum fluctuations (5cr = during inflation, which 

turn into density perturbations after the curvaton dominates the Universe. Unless the ini 



25| 



tial position is close to the hilltop, the WMAP normalization of density perturbations 
imposes the following relation, 

2 

~ 2.4 X 10-^ (12) 



inf 



where H\^i is the Hubble parameter during inflation.'^ Thus, for a given inflation scale 
i^inf, the WMAP normalization determines cr,., whose typical value is /i. Here and in 



^ It is possible to generate density perturbations with a red-tilted spectral index and a vanishingly small 
running, if a initially sits in the vicinity of the hilltop, without relying on large-field inflation. The 
hilltop initial condition generically predicts a sizable non-Gaussianity /nl = 0(10). This enhancement 



of the non-Gaussianity is due to non- uniform onset of curvaton oscillations 2lA l23|. The importance 
of non-uniform onset of oscillations in the curvaton scenario was first pointed out in Ref. |22| to our 
knowledge. 

Precisely speaking, it is necessary to take account of the non-uniform onset of oscillations. However it 
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24| that the predicted density perturbations are close to that in the case of quadratic 



was shown in 

potential unless the initial position is close to the hilltop. 
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what follows we assume that the curvaton dominates the Universe at the decay. We shall 
come back to this issue in Sec. HVl 

As we have seen in the previous section, the observed value of the running spectral 
index indicates that the spectral index is blue at the WMAP pivot scale, and it is red 
at the SPT pivot scale. Thus, the curvature of the potential must be initially positive 
and gradually become negative as the curvaton evolves. Furthermore, the deviation of 
from unity should be comparable to its running. 

First let us see that it is impossible with the curvaton potential (|TTi) to generate such 
sizable running. Using ([9]), ( |T0|) and ( fTTi) . we obtain 

^.-l^^^cos(^), (13) 
a^^^sin^f^V (14) 



9ftHt V/i, 

We note that the spectral index changes from red [ug < 1) to blue [ug > 1), but not the 
other way around. Thus, the running can be only positive, which is also manifest from 
Eq. f|T4|) . In addition, the running is generically of order {rig — 1)^, therefore it is too 
small to account for the observed value j7j. This necessitates some modifications to the 
curvaton potential ffTTl) . 

Let us add modulations 5V{a) to the original potential Vo(cr): 

V{a) = V^{a) + 6V{a), (15) 

with 

5V{a) = h.l(^- cos [j^^y^- (16) 
The relative phase 9 is not relevant for the following analysis, and so, we will set = 0. See 
Fig. [H Such two contributions can be originated from the axion anomalous couplings with 



two different gauge fields, and the multip 
to be ubiquitous in the string axiverse 



e contributions to the axion mass are considered 



18|. Alternatively we may consider a complex 



scalar field $ with an approximate global U(l) symmetry broken by non-renormalizable 
operators. For instance, consider a scalar potential, 

|^|2n-2 / (^r 



V{^) = -m||$|2 + LI + 3 + a^-— 3 + h.c. ) (17) 
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where is a cut-off scale and we assume n < m. Then $ develops a vev at ($) ~ 
(m$M"~'^)^/*^"~^) and the U(l) symmetry gets spontaneously broken. The phase of $ is a 
pseudo-Nambu-Goldstone boson and receives two sinusoidal potentials of different height 
and period from the breaking terms in the parenthesis. Such a scalar potential can be 
realized in a supersymmetric theory with a discrete (R) symmetry. 

For convenience let us introduce two mass scales, = Af/fi, where i runs over 1,2. 
We assume the following hierarchy among the model parameters [l7|. 

» mlfl (18) 
mlh » mlh, (19) 
ml ml, (20) 

which imply /i ^ /2- The first two conditions ensure the curvaton dynamics is basically 
determined by Vq for most of the field range except for the vicinity of the extrema of Vq. 
On the other hand, the third condition tells us that the local curvature of the potential 
is dictated by the modulations, 6V.^ In particular, the WMAP normalization condition 
(|T2|) is maintained, while the spectral index and its running can be significantly affected 
by the modulations 6V. 

Let us now derive conditions on the modulations 6V to realize the observed spectral 
index and the running. From ((Tj) and (1201) . the curvature of 6V should satisfy 

m^cos(^|) ±0(0.01) Hi (21) 

in order to realize — 1 = ± (9(0.01) at the SPT and WMAP pivot scales. On the other 
hand, the slow-roll condition reads 



m 



I < O{0.l)Hl,. (22) 



^ Apart from the running spectral index, such a set-up will be useful when we would like to realize a 
red-tilted spectral index in curvaton scenario without relying on large-field inflation which necessitates 
super-Planckian variation of the inflaton. The modulation generates a red-tilted spectrum, while it is 
Vq that determines the curvaton abundance. 
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FIG. 1: The potential V{a) is shown. We set mi = O.OIH^, m2 = 0.44iJ*, and /2 = 1.25 x 
10~^/i, and the potential and the curvaton field are normalized by mf/f and /i, respectively. 
We can see the modulations 5V is subdominant around the inflection point, but it becomes 
significant around the origin. 



The above two conditions imply 

ml ^ O{0.1)Hl,, (23) 
cos (^j-^ ~ ±0(0.1). (24) 

Furthermore, since the spectral index changes its sign between the two pivot scales, the 
curvaton should vary by about c/2 with c = 0.1 ~ 1 during the e-folds AA^pivot — 2.5. 

Let us evaluate the variation of the curvaton Act while the two pivot scales exited the 
horizon during inflation. The equation of motion can be approximated as 



3H,^i& ~ -mf/isin(^-J , (25) 

where the dot denotes the derivative with respect to time. Since we are interested in the 
variation of |Acr| ^ /2 ^ /i, we may approximate the rhs of (l25l) to be a constant. Thus, 



8 



Ao" is related to the increase of the e-folding number AA^pivot as 



Aa ~ -TTTJI^sin AiVpi^t. (26) 

-^-ninf V /l / 

Combined with Ao" = c/2, we obtain 

|||sin(^)AiVp,™, = ch (27) 

with c ^ 0.1 - 1. 

The decay constant /2 determines the period of the modulations. Both the spectral 
index and its running return to its original value as the curvaton varies by 27r/2. (Here 
we have neglected the change of Vq during one period, because /i ^ /2-) As the SPT 
(and WMAP) data implies an almost constant running over the observed CMB scales, the 
curvaton should not vary by more than the half of the period 71/2 while the CMB scales 
exited the horizon. The corresponding e-folding number is roughly ANqmb — In(i'max) — 
8, where ^max ~ 3000 for the SPT data. Using (l26l) . this constraint reads 

^ sin f ^"j AiVcMB < vr/2, (28) 

which can be satisfied for /2 given by (|27|) . 

Let us write the spectral index and the running in terms of the curvaton potential. 
The spectral index is 

and its running is given by 

Thus we can see that both and a oscillate due to the modulations as the curvaoton 
evolves. In particular, it moves clockwise along an oval in the (n^, a)-plane, as shown in 
Fig. |2j The oscillation amplitudes An^ and Aa are related as 

^ TTT^^n., (31) 
■^^^ pivot 
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where we have used Eq. ([2] 

Similarly, we can easily show that the oscillation amplitude of the running spectral 
index is generally related to that of [us — 1) as 

Aa ^ (32) 

where AA^ denotes the e-folds during one period of the modulations SV. Schematically 
we may express this result as 



27r 

a ~ in^ — 1 



(33) 



where a and Ug — l are understood to represent their typical values. The observed spectral 
index and the running can be explained for AA^ = 20 ~ 30. (The condition (128|1 requires 
AA^/2 > AA^cMB ~ 8.) 

We note that there is an upper bound on AA^ as 



AN 







sin 


4:71 


COS 


(t 


\ns - 


sin 


if 



-1 



< r-^ 7^- (34) 



Therefore, for e.g. Ug — 0.97, AA^ is smaller than a few hundred. In other words, a is 
always larger than of order (n^ — 1)^ in our framework. 

Lastly let us present the numerical results. We set mi = O.Olif*, m2 = 0.44_f/'*, and 
/2 = 1.25 X lO^^/i. The WMAP normalization (fT2l) will determine the relation between /i 
and H^:, and so, all the mass scales are determined if one fixes /i (or H^,). The following 
numerical results do not depend on specific values of /i. The only assumption is that 
the curvaton dominates the Universe before the decay. This condition depends on the 
thermal history of the Universe, especially the reheating temperature, as well as the decay 
rate and the oscillation amplitude of the curvaton. We will discuss this issue in the next 
section. 

In Fig.lll we show the potential V^a) = Vo{a)+6V{a). We can see that the modulations 
do not change much the slope of the potential around the inflection point, while they are 
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significant around the origin. The evolution of the curvaton after the commencement of 
oscillations are considered to be more complicated than without modulations, which will 
be discussed in SeclIVi 

The evolution of Ug and a is shown in Fig. [2l As the curvaton evolves, the predicted 
(n<j, a) moves clockwise along an oval as indicated by the arrow, and it rotates twice 
during about 50 e-folds, i.e., AA^ ^ 25. The solid line represents the numerical result, 
while the the dashed line represents the analytic solution fl29|) and fl30l) . We also show 
the SPT and WMAP pivot scales by the star and the triangle, respectively. The CMB 
scales {C. < 3000) are shown as the thick line. The position of a* is set to be around the 
inflection point of Vq, but the result is not sensitive to the value of cos(cr*//i) unless it 
is extremely close to the extrema of Vq. (Note that the value of cos((7*//2) determines 
the position of the CMB scales along the orbit.) The reason why the actual orbit is not 
a circle is that our analysis relies on the approximation ( !T8|) - ( 120|) . In particular, 5V' 
partially cancels (enhances) when the a is negative (positive). This reduces the change 
of the running if it is negative. We can see from Fig H] that the orbit is slightly flattened 
at the bottom and the running is more or less constant over the CMB scales. We can see 
that the observed spectral index and the running can be indeed realized in our model. 



IV. DISCUSSION AND CONCLUSIONS 



Let us here discuss the evolution of the curvaton after inflation. When the Hubble pa- 
rameter H becomes comparable to mi, the curvaton starts to oscillate with an amplitude 
of order fi. The modulations 5V does not affect this timing even though the curvature 
of the potential is determined by 5V . This was also confirmed by numerical calculations. 
As the Universe expands, the oscillation energy evolves as non-relativistic matter. The 
situation changes when the oscillation energy becomes comparable to the height of the 
modulations 5V ~ '^i/l- The oscillation amplitude at the time is about {1^2/1711) f 2, 
which is much larger than /2, and so, there are many local minima within the oscillation 
amplitude. Then the curvaton will be trapped by one of those minima, and the oscillation 
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a 




FIG. 2: The evolution of the spectral index and the running a. The numerical result is 
shown by the solid line, while the dashed line represents the analytic solution ()29p and (j30p . As 
the curvaton evolves, the predicted {jis^a) moves clockwise along an oval as indicated by the 
arrow, and it rotates twice during about 50 e-folds. We also show the SPT and WMAP pivot 
scales by the star and the triangle, respectively. The CMB scales {I < 3000) are shown as the 
thick line. 

(angular) frequency suddenly changes from m\ to m2. 

One may expect that domain walls are formed during this transition, because the 
curvaton will be trapped by different vacuum at different places^ Actually however the 
domain wall formation is hindered by the bias, i.e., the difference in the energy density 
among those minima, which arises from Vq. To see this, let us first assume the inflaton- 
matter dominated Universe. Then the transition takes place when ~ ^2(/2//i)- The 

^ We expect that the curvaton fluctuations with a characteristic frequency are produced through para- 
metric resonance while the curvaton oscillates in a potential with modulations. Then the curvaton 
oscillations may be spatially inhomogeneous at subhorizon scales. 



12 



bias e between the adjacent vacua will be minimized in the vicinity of the origin as 

e > mlfl (35) 

If the domain walls are formed, the tension is given by ~ 'n^2/|- Assuming that the scaling 
law is achieved soon after the formation, the domain walls annihilate effectively when the 
force of tension acting on the domain walls becomes comparable to the pressure arising 
from the bias, i.e., i^ann ~ mi{mi/m2). However we can see that 

Hf > i/ann =^ mlh > mffu (36) 

which contradicts with ( lT9l) . This means that the domain walls are not formed from 
the beginning because of the large bias in the infiaton-matter domination era. On the 
other hand, if the reheating is completed at the commencement of oscillations, the above 
condition is modified as 

Hf > ifann =^ mlfj > miff, (37) 

which may be satisfied depending on the parameters. For the parameters adopted in the 
numerical calculation, this condition is marginally satisfied, and so, the domain walls will 
annihilate quickly even if they formed. The domain walls may play an important role for 
another choice of parameters, especially in a more general context of the moduli problem 
in the string axiverse. We leave this for future work. In any case the domain walls are not 
likely formed, or even if they are formed, they will soon annihilate in the present case. 

After the reheating of the infiaton, the fraction of the curvaton energy density will in- 
creases until it decays. In order for the curvaton mechanism to work, the curvaton energy 
should be sizable at the curvaton decay. The fraction of the curvaton energy should be 
greater than ~ 0.01 in order not to generate too large non-Gaussianity. This often requires 
a large inflation scale as well as high reheating temperature. Let us comment on how to re- 
alize the curvaton domination. For a fixed inflation scale, the decay constant /i is more or 
less fixed unless we adopt the hilltop initial condition. A smaller mi would delay the com- 
mencement of oscillations, but in the infiaton-matter domination era, it does not increase 
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the energy fraction of the curvaton7 Thus, it may be necessary to reduce the curvaton 
decay rate To-. If the curvaton has an anomalous couphng to massless gauge bosons with 
a strength of 1//?,, the decay rate is considered to be ~ {Nga'^ /256'7i'^) ml/ f^, where a 
denotes the gauge couphng constant, and A''^ denotes the number of decay modes. The 
decay rate wiU be suppressed if the gauge couphng is sufficiently small. Alternatively it is 
possible that the curvaton decays only into a pair of light fermions through the following 
interaction, 

C D m^e'''/^'^p^p, (38) 

where such coupling arises if ipL{R) is charged under the shift symmetry. Then the decay 
rate is suppressed by the light fermion mass as ~ 771^1712/ ff , which will help the 
curvaton to dominate the Universe. 

Here we comment on the difference from Ref. 17|]. First, we have considered a curva- 



ton model, while the large-field infiation was assumed in [17|. Therefore, the tensor mode 
tends to be smaller in the present case, although the present scenario can be straightfor- 
wardly applied to large-field infiation, which produces the tensor mode within the reach of 
future experiments. A relation between a and like f p3|) holds in the case of large-field 
infiation with modulations, if the correction due to the time evolution of H is taken into 
account. On the other hand, it is possible to generate large non-Gaussianity in our model 
if the fraction of the curvaton energy is of (9(0.01 — 0.1) at the decay. One important 
difference is the ease of model building; for the present scenario to work, we need only 
two sinusoidal potential of different height and period, which can be easily realized in a 
usual field theory, and no super-Planckian variation of the scalar field is required. 

So far we have focused on the generation of density perturbations with a large running 
spectral index. In fact we can use the model (fT5|) to generate a red-tilted spectrum 
Us ~ 0.97 with a negligible running, without relying on large-field infiation.^ If there is no 
structure such as modulations, the red-tilted spectrum implies a sizable (tachyonic) mass 
of the curvaton. Thus, the curvaton will start to oscillate soon after infiation ends unless 



^ The fraction of the curvaton energy density increases in the radiation domination for a smaller mi. 
^ Note that the running a is still larger than of order {ug — 1)^, because of ([M]). 
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it initially sits in the vicinity of the hilltop 23|. However, this conclusion can be avoided if 
we add modulations to the curvaton potential. The red-tilted spectrum 1 — Ug = (9(0.01) 
with a negligible running can be realized if we consider large AA^ (i.e., large /2) with 
m\/Hl = (9(0.01). This is because the curvaton dynamics (the onset of oscillations as 
well as the energy density) is determined by Vq, while it is the modulations that affect 
the spectral index and its running. 

In this letter we have shown that the sizable running spectral index can be explained 
in the axionic curvaton model in which the curvaton potential receives two sinusoidal 
contributions of differrent height and period. The spectral index and the running are 
generally related to as Eq. fl33|) . and they are comparable in size for AA^ = 20 ~ 30. Such 
axions with a potential that consists of multiple contributions may be ubiquitous in the 
string axi verse. 
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